One-scale Model for Domain Wall Network Evolution 



o 
o 



> 
(N 

o 
o 

Oh! 



X 



p. p. Avelino,i'2^0 C.J.A.P. Martins,!'^^!] and J.C.R.E. Oliveirai' 

^Centra de Fisica do Porto, Rua do Campo Alegre 687, 4^69-007 Porto, Portugal 
^ Departamento de Fisica da Faculdade de Ciencias da Universidade do Porto, 
Rua do Campo Alegre 687, 4169-007 Porto, Portugal 
^Department of Applied Mathematics and Theoretical Physics, Centre for Mathematical Sciences, 
University of Cambridge, Wilberforce Road, Cambridge CB3 OWA, United Kingdom 

(Dated: 22 July 2005) 

We introduce a new phenomenological one-scale model for the evolution of domain wall networks, 
and test it against high-resolution field theory numerical simulations. We argue that previous numer- 
ical estimates of wall velocities are inaccurate, and suggest a more accurate method of measurement. 
We show that the model provides an adequate approximation to the evolution of key parameters 
characterizing the evolution of the network. We use the model to study possible scaling solutions 
for domain wall networks, and discuss some of their cosmological consequences. 



I. INTRODUCTION 

Topological defects P, are an unavoidable conse- 
quence of cosmological phase transitions. Understanding 
their formation, evolution and cosmological consequences 
is therefore a crucial part of any serious attempt to un- 
derstand the early universe. Most studies to date have 
focused on cosmic strings, which at least in standard sce- 
narios are cosmologically benign, but other elements of 
the defect zoo are also of interest. 

A case in point is that of domain walls. Possibly the 
main reason for their relative neglect is that it was known 
almost from the outset that observational constraints rule 
them out if their symmetry breaking scale is 77 > IMeV 
On the other hand, later on it has been claimed 
that non-standard domain walls can in fact have inter- 
esting cosmological roles, in particular as a realization of 
the so-called 'solid universe' models 0. More recently, 
further motivation for the study of cosmic defects has 
emerged from fundamental theory, namely in the context 
of brane inflation 5, 6J. As far as is presently known, in 
the more realistic such scenarios it is possible to argue 
convincingly, using the Kibble mechanism, that only cos- 
mic strings will form. Nevertheless, one can certainly 
conceive of alternative scenarios where domain walls or 
monopoles could also form 0, . 

Here, after a very brief introduction to domain walls 
(Sect. we extend our recent work on domain wall net- 
works |j9j,0| by deriving (in Sect. IIIII) an analytic model 
for their evolution, in the same spirit of the model of 
Martins and Shellard for cosmic strings PElig. The 
large-scale features of the network are therefore charac- 
terized by a lengthscale (or correlation length) L and 
a microscopically averaged (root-mean-squared) velocity 
V. We will then provide a thorough discussion of the 
possible cosmological scaling solutions of such networks 



(Sect. IIV|) as well as testing the model by comparing 
it with high-resolution field theory Press-Ryden-Spergel- 
type 01 simulations particular, it will be 

shown that the standard (PRS) method of estimating 
wall velocities is inaccurate, and a more reliable method 
will be introduced. Finally (in Sect. IVIj) we will revisit 
some key cosmological consequences of these networks, 
and present our conclusions. 

Since the derivation of the one-scale model for domain 
walls will be done by analogy with that for cosmic strings 
and, as we shall see, their respective evolution equations 
are very similar, we will pursue the analogy further by 
discussing the domain wall and cosmic string contexts 
side by side at various other points along the paper. The 
purpose of this is twofold. First, it shows that for each 
physically allowed regime or scaling solution for the cos- 
mic string case, there will be a qualitatively analogous 
one in the domain wall case. Hence the existing knowl- 
edge about the former case can be helpful in understand- 
ing the latter. But secondly, it will also show that de- 
spite this qualitative similarity there are differences in 
each case at the quantitative level. Indeed, there are 
cases where the different co-dimension leads to substan- 
tial physical differences, and this is why in some cases a 
regime that is cosmologically benign for one defect type 
can be cosmologically disastrous for the other. 

In particular it will be seen that while for cosmic 
strings the linear scaling solution L cx t is an ubiqui- 
tous attractor (both in standard scenarios where strings 
are benign and in non-standard ones where they domi- 
nate the universe — see Ts^), for domain walls the linear 
scaling regime, although possible, is all but irrelevant cos- 
mologically, and universe domination is their long-term 
attractor. Throughout the paper we shall use fundamen- 
tal units, in which c — h — 1. 



II. DOMAIN WALL EVOLUTION 

•Electronic address: )ppavelin@t'c.up.pt| 

tElectronic address: C. J.A.P.Martin s@damtp.cam.ac.uk| Domain walls arise in models with Spontaneously bro- 

tElectronic address: jeolivei@fc. up.ptl ken discrete Symmetries A simple example is that 



2 



of a scalar field </> with the Lagrangian 

^ = ^'/',a</''" - , (1) 

where the potential V((j)) has a discrete set of degenerate 
minima, say for example 

n0) = ^o(|-i)'. (2) 

By varying the action 

5 = J dtj d^xyf^C (3) 
with respect to we obtain the field equation of motion 

where V is the Laplacian in physical coordinates and 
H = {da/dt)/a is the Hubble parameter. 

In many cosmological contexts of interest, one can 
neglect the domain wall thickness when compared to 
its other dimensions, and thus treat the wall as an in- 
finitely thin surface. With this assumption, its space- 
time history can be represented by a 3D worldshect 
— x^^iC^"-)^ a = 0, 1, 2. A new action can then be eas- 
ily derived, see for example 0|- In the vicinity of the 
worldsheet a convenient coordinate choice is the normal 
distance from the surface. Noticing that in the thin wall 
limit all fields in the Lagrangian should depend only on 
this normal coordinate, and integrating out this depen- 
dence, one finds 

S = -ajd\^, (5) 

where 

lab = gt,vX^^aX'',b (6) 

is the worldsheet metric, with the obvious definition 7 — 
det{iab), and a is the mass per unit area of the wall. 
Notice that this action is proportional to the 3- volume of 
the wall's worldsheet, and hence is clearly the analogue 
of the Goto-Nambu action for strings. Corrections to 
the action due to the finite width of the wall have been 
discussed in T^. 

The equations of motion for a domain wall can then be 
derived by the well-known process of varying the action, 
yielding in this case 

i-"^da (V77'^H) + r(;,7"''^ra-<6 = . (7) 

Since the wall action is invariant under worldsheet re- 
parametrisation, we are free to impose three arbitrary 
gauge conditions. However, no such conditions have been 
found that will lead to equations that can be readily 
solved, either analytically or numerically, as in the case of 



cosmic strings — see [l7| for an approximate analysis. As 
a consequence, much less is known about the dynamics 
of domain walls than about the dynamics of strings. In 
passing, however, we mention a curious point of contact 
between the two cases: in fiat (Minkowski) spacetime, 
any planar string solution can be trivially turned into a 
domain wall solution by simply translating in the direc- 
tion perpendicular to the plane — see |^ for a discussion. 

III. THE ONE-SCALE MODEL 

We now concentrate on analytic modelling, and start 
by presenting a simple phenomenological derivation of a 
one-scale model for defect evolution. We will see that this 
easily reproduces the broad features of one-scale mod- 
els for cosmic strings Q m III El (although strictly 
speaking it is not self-sufficient in this regard), but has 
the advantage of allowing a straightforward extension to 
the case of domain walls. Other approaches to the an- 
alytic modelling of cosmic strings and domain walls can 
be found in 18,, ,19]. 

Let us start by considering a network of non- 
interacting straight infinite strings, oriented along a fixed 
direction and all with the same value of the velocity 
in a fiat FRW universe. Since the momentum per unit 
comoving length is proportional to we have W7 oc a^^ 
or equivalently 

^ + 2i7(l-i;2)^ = 0, (8) 

where t is the physical time, a is the scale factor, H = d/a 
is the Hubble parameter and 7 = (1 — v'^)~^^'^ . On the 
other hand, since the average number of strings in a fixed 
comoving volume should be conserved, one has p oc 7a^^ 
or equivalently 

^ + 2i/(l + i;2)p^0, (9) 

where p is the average energy density in cosmic strings 
and we have also used Eqn. ||SJ| to obtain 0. Note that 
although we are assuming the strings to be infinite and 
straight we may define a characteristic length scale, 

= ^, (10) 
P 

which is directly related to the average distance between 
adjacent strings. Here we have also defined a string mass 
per unit length, ^. 

Of course, this case with non-intersecting straight infi- 
nite strings all aligned along a fixed direction and having 
the same value of the velocity is completely unrealistic. 
In practical situations the special length scale L, defined 
above, will also be approximately equal to the curvature 
scale of the strings, and the strings will have a non-zero 
probability of crossing and interacting with each other. 
Also the value of the velocity will vary along the strings 
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but in the absence of interactions Eqn. (jS)) would re- 
main valid, with v being taken as the RMS velocity of 
the strings. 

Hence, we need to add further terms to Eqns. ^ and 
© in the context of realistic models. Let us consider the 
latter first. The probability of a string element of size L 
encountering one other segment of the same size within a 
time dt is proportional to vdt/L . We may therefore add a 
new term to the right hand side of © , which will account 
for the energy lost from the long string network due to 
loop production, so that the equation now becomes 



or, writing it in terms of the length scale L, 



dL 
) 

' dt 



2(1 



)HL 



(11) 



(12) 



This is of course under the assumption that strings do in- 
tercommute when they interact — see for a discussion 
of alternative scenarios. On the other hand, in the above 
we assumed straight infinite strings. Realistic strings will 
of course be curved, and their curvature will be responsi- 
ble for an acceleration term which also needs to be taken 
into account. The velocity equation is thus corrected to 



dv 



(1 



2Hv\ 



(13) 



Notice that (|12I13|) are the evolution equations of the 
simplest version of the VOS model [ij, the only 

things missing are the fact that the curvature correction 
term k should be velocity-dependent jl3j , and additional 
terms accounting for effects like friction due to particle 
scattering, spatial curvature and so on. It is worth stress- 
ing at this point that one of the basic assumptions of a 
one-scale model is, not surprisingly, that there is a sin- 
gle large relevant lengthscale in the problem. In addition 
to the characteristic lengthscale L defined above (which 
is essentially a parametrization of the energy density in 
strings) one can define a correlation length ^ and a curva- 
ture radius R, for example, and one can certainly envis- 
age them being different — in fact this can be confirmed 
numerically J23 , and other approaches to modelling do 
allow for it [ig. However, from the point of view of one- 
scale modelling, one is assuming that L = ^ = i?. 

Despite the simplicity of the argument, the above 
derivation shows that this reasoning is robust, which is 
useful since we can easily apply it to the case of domain 
wall networks. Again we start by considering a network 
of non-interacting infinite planar domain walls oriented 
along a fixed direction and all with the same value of 
the velocity v in a flat FRW universe. The momentum 
per unit comoving area is proportional to so that we 
have U7 oc or equivalently 



dv 
'dt 



+ '3H{l-v^)v = 0. 



(14) 



As before we will assume that the average number of 
domain walls in a fixed comoving should be conserved so 
that p oc 7a^^ or equivalently 



(15) 



where p is now the average energy density in domain 
walls and we have again used Eqn. 114|) to obtain H15|l . 
Similarly, although we are assuming the domain walls 
to be infinite and planar we may define a characteristic 
length scale, 



P 



(16) 



which is directly related to the average distance between 
adjacent walls measured in the frame comoving with the 
expansion of the universe, and a is now the domain wall 
mass per unit area. Finally, making analogous modifica- 
tions to allow for energy losses and acceleration due to 
the wall curvature, we find 



dL 



(1 + 'iv^)HL + 



CwV 



and 



|^(l-.^)(^_3i.. 



(17) 



(18) 



the latter equation has also been previously obtained, 
using a different approximation, in [l7|. 

Note that there is a difference between the two cases 
when it comes to energy losses by intercommuting. String 
loops are dynamically very important due to the ex- 
istence of a range of trajectories that are not self- 
intersecting, and hence are long-lived. No such solu- 
tions are known, or believed to exist, for domain walls, 
so whenever closed walls (also called 'vacuum bags') are 
produced they will decay very quickly — a point already 
made by . 

These therefore provide a phenomenological model for 
domain wall evolution. As in the case of strings, we 
may hope that the energy loss efficiency Cw may be 
constant and possibly independent of the cosmological 
epoch, but by the same token we expect the curvature 
parameter to be a velocity-dependent function. In 
the case of cosmic strin gs, one can use a combination 
of field theory [23, [23 and Goto-Nambu numeri- 
cal simulations |23, |2J| to provide a successful calibra- 
tion, leading to the VOS model [23, [l^- For domain 
wa lls, however, only field theory simulations are available 
IlEflllllllllllMlllIll (though a simulation in the 
thin wall approximation has been attempted in JXl), so 
no similarly accurate calibration will be possible. Still we 
will show that as far as one can ascertain, the model does 
provide an adequate description of our numerical simu- 
lations. Before this however, we will explore the possible 
scaling solutions of the model. 
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IV. SCALING SOLUTIONS 



We now use the evolution equations (|17I18() of our phe- 
nomenological model to discuss the cosmological evolu- 
tion of domain wall networks, in particular discussing all 
relevant scaling solutions. The derivation of these scaling 
laws is mostly analogous to what has been done for cos- 
mic string networks in the case of the VOS model. Hence 
we will present each solution in turn and discuss the phys- 
ical mechanisms behind them, but will not in general 
present detailed derivations. We refer the reader to the 
original papers on cosmic string networks ^> ^> ^ 
for more detailed discussions. Most of the derivations are 
fairly straightforward, even if somewhat tedious. 

Let us start by neglecting the effect of the energy den- 
sity in the domain walls on the background (specifically, 
on the Friedmann equations). As we shall shortly see this 
is not a good approximation, since the wall network will 
generally end up dominating the energy density of the 
universe. However it is this scenario that is effectively 
considered, for example, when one performs numerical 
simulations of domain wall networks. 

In this case it is easy to see that, just as for cosmic 
string networks, the attractor solution to the evolution 
equations (|17I18|I corresponds to a linear scaling solution 



L^et. 



V — const . 



(19) 



Assuming that the scale factor behaves as a cx the 
detailed form of the above linear scaling constants is 



3a(l — a) 



1 — a 
3a ku 



(20) 



(21) 



As in the case of cosmic strings ]^ , an energy loss mech- 
anism (that is, a non-zero c^) may not be needed in order 
to have linear scaling. In fact, by considering the Q 
limit one finds that for 
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(22) 



a linear scaling solution is always possible. Hence in this 
case a linear scaling solution may exist in both mat- 
ter and radiation eras (in the case of cosmic strings 
this is only guaranteed to be the case in the matter 
era.) This means that having non-standard (that is non- 
intercommuting) domain walls is by no means sufficient 
to ensure a frustrated wall network. 

On the other hand, if a < 1/4 then an energy loss 
mechanism is necessary to have linear scaling. Note, how- 
ever, that the linear scaling solutions are physically very 
different for cosmic strings and domain walls. In the case 
of cosmic strings, in the linear scaling phase the string 
density is a constant fraction of the background density, 
whereas in the case of domain walls we have 



CX t . 



Pb 



(23) 



so the wall density grows relative to the background den- 
sity, and will eventually become dominant. This happens 
at a time 



U ^ (GaY 



(24) 



Since the domain wall mass per unit area is related to 
the energy scale of the phase transition, a ^ r]^, we can 
also write out a given epoch as 

tpi \mpi 



(25) 



hence walls that would become dominant around today 
would have been formed at a phase transition with an 
energy scale 



r]o ~ lOOAfeF; 



(26) 



notice that this is two orders of magnitude larger than 
the standard Zcl'dovich-Kobzarev-Okun bound j^]. It 
will be seen from the discussion that follows that net- 
works that are much heavier would have become domi- 
nant well before having reached the linear scaling regime, 
whereas networks that are much lighter would not yet 
have reached the linear regime by today. Hence the range 
of cosmological scenarios where the linear scaling solution 
is of interest is quite limited. 

There is, moreover, an effect which we have neglected 
thus far. At early times, in addition to the damping 
caused by the Hubble expansion, there is a further damp- 
ing term coming from friction due to particle scatter- 
ing off the domain walls. Phenomenologically, it can be 
shown 2] that its effect can be adequately described by 
a frictional force per unit area 



(27) 



where v is the string velocity. In the above we have de- 
fined a friction length scale 



CX a 



(28) 



where T is the temperature of the background and is 
the number of light particles changing their mass across 
the walls 1]. If the self-coupling of the domain wall field 
and its couplings to the other fields are not very small, 
then this expression holds at all times after the wall for- 
mation. In other cases the behaviour might be slightly 
different very close to the phase transition 0, but since 
we are mostly interested in the behaviour of wall net- 
works at recent times we shall neglect this subtlety. (Also 
can be effectively zero at low temperature if the wall 
is described by a single field.) 

Just like in the case of cosmic strings 0, 0, is 
then easy to modify the evolution equations of our one- 
scale model to account for this extra friction term. They 
become 



dL 



L 2 

HL + —V + CyjV 



(29) 
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where we have defined a damping length scale 

which includes both the effects of Hubble damping and 
particle scattering. If no particles scatter off the walls 
{Nu, — 0) then the friction length scale is infinite, and the 
only damping term comes from Hubble damping. Note 
that since £/ oc a^, the friction term will be dominant 
at early times, while the Hubble term will dominate at 
late times, so the late-time linear scaling solution is un- 
changed. However, notice that the timescale when Hub- 
ble damping dominates over friction (which is also the 
timescale for the walls to become rclativistic) is again t^, 
given by Eqn (|24|l . Thus we see that domain wall net- 
works will dominate the energy density of the universe 
even without ever becoming relativistic or reaching the 
linear scaling regime. 

Just as in the case of cosmic strings 0, 0| , there will 
be two possible scaling solutions (which are necessarily 
transient) during the friction-dominated epoch. Also as 
in the case of strings these solutions will exist re- 
gardless of whether or not the walls interact with each 
other (that, is, whether Cw is non-zero or vanishes). If the 
defect-forming phase transition is such as to produce a 
low-density network there will be an initial period where 
the network will be conformally stretched. The scaling 
laws will therefore be 



oc a , 



Vs oc 



(32) 



(33) 



Notice that for domain walls this gives w oc a'^, whereas 
for cosmic strings we would have w oc a^. Going back to 
the domain walls we respectively have, in the radiation 
and matter-dominated epochs 



Lr (X t 



1/2 



3/2 



and 



Lrn OC t 



2/3 



Vr oc t 



Vm Oct 



(34) 



(35) 



We emphasize that although the network is being 
stretched as the scale factor, and is non-relativistic, the 
velocities are increasing rather fast, due to the effect of 
the domain wall curvature. This indicates that even in 
the absence of other mechanisms this regime would only 
be a transient. The only situation where such a stretch- 
ing regime could persist would be during an inflationary 
phase, but in that context the much faster expansion is 
enough to counter the wall velocities, and in fact make 
them decrease. It's indeed easy to see that in the case of 



an exponential expansion the solution of H29I30|I has the 
form 



-'inf 



Vinf 



(36) 



(37) 



Following the conformal stretching regime, or perhaps 
right after the formation of the network if it is formed 
with high enough density, there is a further transient 
scaling regime. The inevitability of such a regime for 
both strings and walls was first argued for, using simple 
physical arguments, in T|, so we shall call this the Kibble 
regime. In the context of velocity-dependent models the 
existence of the Kibble regime can be rigorously derived. 
The scaling solution has precisely the same form for both 
types of defects 



Lk oc 



1/2 



Vk OC {ifHf'^ 



(38) 



(39) 



although of course the friction lengthscale will not have 
the same form in the two cases. For domain walls we 
respectively have, in the radiation and matter-dominated 
epochs 



Lr oc t^/'^ , Vr oc t^^^ 



and 



Lm, oc t 



11/6 



Vm. OC t 



5/6 



(40) 



(41) 



Notice the differences relative to the stretching regime. 
Here the correlation lengths grow much faster, while the 
velocities grow relatively more slowly. Physically, the 
difference between the two regimes is one of interactions 
and energy losses. In the stretching regime the walls are 
typically quite far apart, so there is very little interac- 
tion between them — typically less than one per Hubble 
volume per Hubble time. In the Kibble regime, on the 
other hand, the walls are so close together that there 
is a very large number of interactions — in fact there are 
more than in the case of the linear scaling regime. This 
enhanced energy loss makes the correlation length grow 
quite fast. The wall velocities are still non-relativistic and 
growing, but because regions of the network with higher 
velocity than average have a larger interaction probabil- 
ity than slower regions (theupon leaving the network) the 
enhanced energy loss is also responsible for making the 
velocities grow more slowly than in the stretching case. 
Still the Kibble scaling is also a transient, which in the 
absence of other mechanisms will necessarily end when 
the network becomes relativistic. 

In passing, it is again interesting to compare these scal- 
ing solutions to the ones for cosmic strings, as derived in 
Il3l. We notice that wall velocities become relativistic 
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faster than those of strings, and also that wall densities 
grow faster relative to the background. In fact in both 
the stretching and Kibble regimes the wall density grows 
relative to the background. On the other hand, the string 
density grows relative to the background in the stretching 
regime, but decreases relative to the background in the 
Kibble regime. Thus the friction-dominated epoch will 
last comparatively less for domain walls than for strings. 

Even allowing for friction, linear scaling would be an 
attractor of the above equations if one neglected the ef- 
fect of the wall density on the expansion of the universe. 
However, we have seen that in every scaling regime con- 
sidered the wall density grows relative to the background, 
so that a wall density term 

P» = ^ (42) 

must be included in the Einstein equations. This changes 
the situation for it is easy to see that the domain wall 
network will eventually dominate the energy density of 
the universe (unless some mechanism like a subsequent 
phase transition were to make it decay and disappear). 

Thus we again see that linear scaling is of little practi- 
cal importance, since it is never reached for any cosmolog- 
ically realistic network. Heavy domain walls will quickly 
dominate the energy density of the universe, thus chang- 
ing the behaviour of the Friedmann equation, typically 
before linear scaling is reached — in any case any networks 
where it could have been reached already must be suffi- 
ciently heavy to be observationally ruled out. Light walls, 
on the other hand, won't yet have reached that solution: 
their dynamics will still be friction-dominated today if 
they are to be cosmologically viable. 

Since a domain wall network will eventually dominate 
the energy density of the universe it is important to study 
the dynamics of the universe in this regime. The expecta- 
tion 13 is that the domain wall network will again become 
frozen in comoving coordinates with 

L(xa (43) 

so that the scale factor should now grow as 

a(xt^ . (44) 

In this case the average distance between the walls also 
grows as and rapidly becomes greater than the hori- 
zon. This will happen at a time that is again given by 
ti, in Eqn. (|24|l . An inertial observer will see domain 
walls moving away towards the horizon, and as walls fade 
away the spacetime around the observer will asymptot- 
ically approach Minkowski space. Notice that this solu- 
tion does not depend on Cw — it is valid whether or not 
the domain walls interact. In the case of cosmic strings 
it has been shown 0, 0| that the asymptotic solution 
is always Ls oc t (again, whether or not the strings in- 
tercommute). The onset of domain wall domination will 
be studied in more detail (both analytically and numer- 
ically) in future work. 



V. NUMERICAL SIMULATIONS 

We now proceed to compare the predictions of our ana- 
lytic phenomenological model, whose evolution equations 
are (|17I18|I . with the results of numerical simulations. 
We shall use our own set of high-resolution field theory 
numerical simulations, done on the COSMOS supercom- 
puter and using the PRS algorithm > which were first 
discussed in As we said above, there is no ana- 

logue of the Goto-Nambu ID simulations for the case of 
domain walls, though an attempt in that direction has 
been performed (with moderate success) in [T^ . 

The fact that linear scaling is of almost no relevance 
in realistic cosmological scenarios involving domain walls 
counts against producing a well-calibrated model here, 
for the scenarios that one can easily simulate numerically 
are not the most interesting ones in practice, and vice- 
versa. A further difficulty with the calibration of the 
code using simulations with the PRS algorithm is 
that this modifies the wall thickness. This has the side 
effect of erasing small-scale structures on the walls, and 
hence also destroying small closed walls when the wall 
thickness becomes comparable to their size. 

A. Measuring wall velocities 

A crucial point, which has been somewhat neglected 
in previous analyses, is that estimating the velocities 
of domain walls is notoriously difficult in field theory 
simulations — see for a detailed discussion of the prob- 
lems involved in the case of cosmic strings. Earlier field 
theory simulations, using the PRS algorithm jl^J typi- 
cally find 

Vprs ~ 0.4 , (45) 

no significant difference being found (when allowance is 
made for the magnitude of the numerical errors involved) 
in the values for 2D and 3D simulations. On the other 
hand, Kawano |l7j . which as we have said uses a thin wall 
approximation (akin to the Goto-Nambu simulations for 
strings), and further simplifies the problem by consider- 
ing only 2D simulations, typically finds smaller velocities 

vkaw ~ 0.25 - 0.30 . (46) 

We note that also in the case of cosmic strings it has been 
observed that velocities measured from field theory simu- 
lations tend to be somewhat higher — again see [23| for a 
detailed discussion. Both types of simulations find some 
evidence for an approach to linear scaling, both in the 
radiation and matter eras. A more detailed discussing of 
the approach to linear scaling in domain wall field the- 
ory simulations can be found in fl3| . These findings 
prompt us to study the issue of velocity estimations in 
more detail. 

The frictionless evolution of a domain wall network in 
a flat homogeneous and isotropic Friedmann-Robertson- 
Walker (FRW) universe is given by Eqn. with the 
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potential given by Eqn. If we neglet the damping 

term then the solution for a plane wall with velocity v is 
given by 



tanh 



7r7o 



(z - zo - vt) 



(47) 



1 /2 

where z is a physical coordinate and uj = 7r(^o/(2Vb) 
is the constant wall thickness in physical coordinates. It 
is straightforward to show that the ratio between the ki- 
netic part of the wall energy density and the total energy 
density p(z) is independent of the z coordinate and is 
equal to 



p{z) 2^lV{4>) ' 



where 



|V</ 



(48) 



(49) 



The domain walls in realistic network simulations will 
obviously not be planar and the velocities will vary along 
the walls. However, an estimate of the microscopic rms 
velocity in field theory simulations can be made using 
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(50) 



where the sum is limited to the points on the grid which 
intersect a domain wall. We define the wall region to be 
the region for which — e^o < 4> < e4>o but of course there 
is some ambiguity on the choice of e. However, taking 
into account the result obtained for the planar wall one 
expects that the results will converge for small enough e 
(note that we cannot realisticaly make e arbitrarily small, 
since in that case we would have very poor statistics). 
This is indeed what we find. Our approximation almost 
completely eliminates the radiated energy from the walls 
which otherwise would contaminate the estimate of the 
velocities. This is an important advantage over previ- 
ous estimations of v where the energy radiated from the 
walls was not eliminated in the velocity estimations. We 
therefore expect that these previous estimates have over- 
estimated the wall velocities. 

Following we have modified the equations of mo- 
tion in such a way that the walls' co-moving thickness 
is fixed in co-moving coordinates in order to be able to 
resolve the domain walls through the network evolution. 
Ignoring the damping term due to the expansion of the 
universe the modified equations become 



drf 



V 



dy_ 



(51) 



where is the Laplacian in comoving coordinates. The 
planar wall solution can now be written as 



bo tanh 



7r7 



(52) 




FIG. 1: Comparing the standard (PRS) and our (new) ve- 
locity estimation methods, for 2D (solid) and 3D (dashed) 
radiation era runs. Each curve corresponds to an average of 
10 simulations. Clearly the PRS method overestimates veloc- 
ities by about 20 percent. 



Notice that here w is a fixed comoving thickness so that 
the physical wall thickness decreases with time. Con- 
sequently, with this approximations the estimation of 
{v^^'^) in realistic domain wall network simulations must 
now be made as 



21^(0) 



rather than 



(53) 



(54) 



It is easy to show from Eqn. Q that for a planar wall 
solution 



dt a 



0, 



(55) 



so that (0) oc and 7U oc . In order for the momen- 
tum conservation law of the wall evolution in an expand- 
ing universe to be maintained one must add a damping 
term to Eqn. H51|) so that now {dcfi/drj) oc and we 
have 



d^(t) 3 da d(j) 
drf a drj di] 



dV 



(56) 



The standard (PRS) method of estimating velocities 
can be readily compared with ours. We have thus es- 
timated the velocities by both methods in series of 2D 
and 3D, radiation and matter era runs. Velocities were 
measured at all timesteps, and then sorted in increasing 
order, so that the two velocities can be ploted directly 
against each other without reference to time. A stan- 
dard cloud-in-cell algorithm was used for data smooth- 
ing. The results of this analysis are shown in Fig. ^and 
121 respectively for the radiation and matter eras. 
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FIG. 2: Same as Fig. Qfor matter era simulations. 

We can see that the standard PRS method of veloc- 
ity estimation typically overestimates velocities by some- 
thing like 20 percent. Note that in either case the result 
is the same for 2D and 3D simulations (when allowance 
is made for numerical error bars, which were not plot- 
ted), and there is also very little difference between the 
radiation and matter cases. Indeed, if one fits a linear 
function to the above data, we find for the radiation era 

Vrad = 0.79 Vprs + 0.03, (57) 




FIG. 3: Comparing our analytic model with field theory nu- 
merical simulations in the matter (left-side panels) and radia- 
tion (right-side panels) epochs. Top panels show the evolution 
of the correlation length (in fact L/t) whereas bottom panels 
show the RMS velocity. In both cases the longer solid curves 
correspond to 2D simulations, while the shorter solid curves 
correspond to 3D simulations. Each such curve is in fact an 
averaged curve for the outcome of 100 such simulations. The 
corresponding dashed curves are the outcome of our analytic 
model, for the initial conditions relevant to each case. For 
all such curves we have used Cm ~ 0.5 and fcm — 0.9, which 
provides the best fit. 



while in the matter era 

Vmat^0.83Vprs+0.00. (58) 

Notice that these fitting functions should not be used out- 
side the specified range. In particular, one expects that 
there will be deviations from this simple linear behaviour 
for very small (non-relativistic) velocities. 



B. Testing the analytic model 

Mindful of these caveats, we have compared our an- 
alytic model with the numerical simulations of [lol |. 
We took four series of high-resolution field theory simu- 
lations, in 2D and 3D and for the radiation and matter 
dominated epochs. Each such series is composed of 100 
different runs. For each series we have averaged the re- 
sult of the 100 runs, in particular calculating averaged 
correlation lengths and RMS velocities that can be read- 
ily compared with the predictions of the model. As we 
have pointed out, although we expect the parameter 
to have some dependence on velocity, it is not easy to 
determine it. Hence, as a first approximation, we shall 
start by assuming that it is a constant, just like c^,. In 
these circumstances, we find that the best fit to the sim- 
ulations is the one shown in Fig. 13 which corresponds 



to the parameters 

c^, ~ 0.5 , ~ 0.9 . (59) 

Several comments are in order here. The first one is 
that the fit is quite good (especially for the velocities), 
given that as we said we know that a constant k^ is 
only an approximation. Importantly, the same parame- 
ters provide a good fit both in the radiation and in the 
matter epochs. (Looking for separate best fits in the two 
epochs would lead to only slightly different parameters, 
but given the inherent error bars they would effectively be 
indistinguishable.) Interestingly, the fit would be much 
worse if we had used the velocities measured by the stan- 
dard PRS method. 

Obviously the fit is only good at late times. At early 
times, one is still very close to the phase transition and 
fields are relaxing, so the domain wall network is not yet 
well-defined, and obviously the analytic model cannot be 
expected to account for such complicated dynamics. We 
could nevertheless improve the early-time fit if we had al- 
lowed for friction and radiation terms in the model — this 
has been successfully done for field theor y si mulations of 
cosmic string networks, as described in j23- Given the 
simplicity of the present model and the approximations 
being made elsewhere, we think this would be an un- 
necessarily complication at this stage, though it should 
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t t 

FIG. 4: An alternative fit to the model, now using an ansatz 
for kwiv) that has been derived for cosmic strings — see main 
text for further discussion. Panel and line conventions are as 
in Fig. H 



certainly be addressed in the future. 

We emphasize that we are dealing with a phenomeno- 
logical model (in fact, rather more so than in the case 
of cosmic strings), so one should not attach too deep a 
meaning to the parameters one finds. Presumably these 
will change if one has a proper ansatz for kyj. Even 
though we have not attempted to calculate this, we can 
for the sake of the argument see what happens if we as- 
sume that the function k{v) that was obtained for cosmic 
string networks is also applicable for domain walls. This 
has the following form 

k{v) = ^(1 - -^)(i + 2V2.3)^q:|5 ; (60) 

we refer the reader to ^3 fo'^ ^ detailed derivation and 
further discussion. Notice that with this ansatz k is no 
longer a free parameter, and it is always smaller than 
unity. The result of this alternative fit is shown in Fig. 
0] and we now find a best-fit parameter 

c™~1.0. (61) 

As expected, this does not provide a very good fit. 
The energy loss term is now much larger, which has an 
effect on the velocities. Again a single values provides a 
reasonable fit for the correlation lengths in the radiation 
and matter epochs, but the new fit significantly underes- 
timates the velocities. The velocity fit could of course be 
improved, at the expense of a poor fit for the correnation 
lengths. 



VI. CONCLUSIONS 

We have introduced a new phenomenological one-scale 
model for the evolution of domain wall networks, by anal- 
ogy with the stardard analytic model for cosmic string 
evolution. This model has then been tested against high- 
resolution field theory numerical simulations (in radia- 
tion and matter dominated epochs), and a good agree- 
ment has been found. Importantly, we have argued that 
previous numerical estimates of wall velocities are inac- 
curate, and tend to overestimate the wall velocities due 
to the effect of the radiation background. We have quan- 
tified the inherent inaccuracy (which is found to be at 
around the twenty percent level) and provided a more 
accurate method of measurement. Our phenomenological 
model provides an adequate approximation to the evolu- 
tion of key parameters characterizing the evolution of the 
network, but more accurate modelling can lead to even 
better fitting. 

We have also used the analytic model to exhaustively 
study possible scaling solutions for domain wall networks, 
and discuss some of their cosmological consequences. 
Comparison with analogous results for cosmic string net- 
works has led to the identification of a number of similar- 
ities and differences between the two cases. Indeed, for 
each physically allowed regime or scaling solution for the 
cosmic string case, there is be a qualitatively analogous 
one in the domain wall case. However, beyond this qual- 
itative similarity there are differences in each case at the 
quantitative level. The different co-dimension leads to 
substantial physical differences, and this is why in some 
cases a regime that is cosmologically benign for one de- 
fect type can be cosmologically disastrous for the other. 

A case in point is that while for cosmic strings the 
linear scaling solution L oc t is an ubiquitous attractor 
(both in standard scenarios where strings are benign and 
in non-standard ones where they dominate the universe), 
for domain walls the linear scaling regime, although pos- 
sible, is all but irrelevant cosmologically, and universe 
domination (with L oc t^) is their long-term attractor. 

Finally, let us end by noting that in the present work 
we have restricted ourselves to numerical simulations of 
domain walls in universes dominated by radiation and 
matter. These are of course the easiest scenarios to sim- 
ulate, but plainly they are not the most relevant cos- 
mologically. It would therefore be interesting to carry 
out analogous simulations of domain-wall dominated uni- 
verses, not unly to understand how the domination sets 
in but also to test our phenomenological model in this 
context. We will tackle this issue in a future publication. 



VII. ACKNOWLEDGMENTS 

This work was funded by FCT (Portugal), in 
the framework of the POCI2010 program, sup- 
ported by FEDER. Specific funding came from grant 
POCTI/CTE-AST/60808/2004 and from the Ph.D. 



10 



grant SFRH/BD/4568/2001 (J.O.) 

This work was done in the context of the ESF COSLAB 
network. Numerical simulations were performed on 



COSMOS, the Altix3700 owned by the UK Computa- 
tional Cosmology Consortium, supported by SGI, Intel, 
HEFCE and PPARC. 



T. W. B. Kibble, J. Phys. A9, 1387 (1976). 

A. Vilenkin and E. P. S. Shellard, Cosmic Strings and [17 
other Topological Defects (Cambridge University Press, [18 
Cambridge, U.K., 1994). 

Y. B. Zeldovicii, 1. Y. Kobzarev, and L. B. Okun, Zh. [19 
Eksp. Teor. Fiz. 67, 3 (1974). 
M. Bucher and D. N. Spergel, Phys. Rev. D60, 043505 [20 
(1999), astro-ph/9812022. 

N. Jones, H. Stoica, and S. H. H. Tye, JHEP 07, 051 [21 

(2002), hep-th/0203163. 

S. Sarangi and S. H. H. Tye, Phys. Lett. B536, 185 [22 
(2002), hep-th/0204074. 

T. Matsuda, JHEP 10, 042 (2004), hep-ph/0406064. [23 

N. Barnaby, A. Berndsen, J. M. CUne, and H. Stoica, 

JHEP 06, 075 (2005), hep-th/0412095. [24 

J. C. R. E. Oliveira, C. J. A. P. Martins, and 

P. P. Avelino, Piiys. Rev. D71, 083509 (2005), hep- [25 

ph/0410356. 

P. P. Avelino, J. C. R. E. Oliveira, and C. J. A. P. Mar- [26 
tins, Phys. Lett. B610, 1 (2005), hep-th/0503226. 
C. .J. A. P. Martins and E. P. S. Shellard, Phys. Rev. [27 
D53, 575 (1996), hep-ph/9507335. 

C. J. A. P. Martins and E. P. S. Shellard, Phys. Rev. [28 
D54, 2535 (1996), hep-ph/9602271. 

C. J. A. P. Martins and E. P. S. Shellard, Phys. Rev. [29 
D65, 043514 (2002), hep-ph/0003298. 
W. H. Press, B. S. Ryden, and D. N. Spergel, Astrophys. 
J. 347, 590 (1989). [30 
C. J. A. P. Martins, Phys. Rev. D7G, 107302 (2004), 
hep-ph/0410326. 

R. Gregory, D. Haws, and D. Garfinkle, Phys. Rev. D42, 



343 (1990). 

L. Kawano, Phys. Rev. D41, 1013 (1990). 

D. Austin, E. J. Copeland, and T. W. B. Kibble, Phys. 

Rev. D48, 5594 (1993), hep-ph/9307325. 

M. Hindmarsh, Phys. Rev. D68, 043510 (2003), hep- 

ph/0207267. 

J. N. Moore, E. P. S. Shellard, and C. J. A. P. Martins, 
Phys. Rev. D65, 023503 (2002), hep-ph/0107171. 
U.-L. Pen, U. Seljak, and N. Turok, Phys. Rev. Lett. 79, 
1611 (1997), astro-ph/9704165. 

G. Vincent, N. D. Antunes, and M. Hindmarsh, Phys. 

Rev. Lett. 80, 2277 (1998), hep-ph/9708427. 

D. P. Bennett and F. R. Bouchet, Phys. Rev. D41, 2408 

(1990). 

B. Allen and E. P. S. Shellard, Phys. Rev. Lett. 64, 119 
(1990). 

C. J. A. P. Martins, J. N. Moore, and E. P. S. Shellard, 
Phys. Rev. Lett. 92, 251601 (2004), hcp-ph/0310255. 

D. Coulson, Z. Lalak, and B. A. Ovrut, Phys. Rev. D53, 
4237 (1996). 

S. E. Larsson, S. Sarkar, and P. L. White, Phys. Rev. 
D55, 5129 (1997), hep-ph/9608319. 

P. P. Avelino and C. J. A. P. Martins, Phys. Rev. D62, 
103510 (2000), astro-ph/0003231. 

P. P. Avehno, J. P. M. Carvalho, C. J. A. P. Martins, 
and J. C. R. E. Oliveira, Phys. Lett. B515, 148 (2001), 
astro-ph/0004227. 

T. Garagounis and M. Hindmarsh, Phys. Rev. D68, 
103506 (2003), hep-ph/0212359. 



